Photoassociation spectroscopy of ^^Sr: reconstruction of the wave function near the 

last node 



O 

o 



Oh, 



O ■ 

Oh: 



(N 

> 
m 
in 
o 

o 

o 

o 



X 



Masami Yasuda,* Tetsuo Kishimoto, Masao Takamoto, and Hidetoshi Katori 

Engineering Research Institute, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan 

(Dated: February 2, 2008) 

We have investigated photoassociation (PA) spectra of ultracold **Sr atoms near the ^So + 
5s5p ^Pi atomic asymptote. The intensity modulation of the PA hues was used to reconstruct the 
ground-state scattering wave function, whose last nodal point was determined to be rg = 3.78(18) 
nm. The PA lines also determined a precise lifetime of the ^Pi state to be 5.263(4) ns. 

PACS numbers: 34.50.Rk, 32.70.Cs, 32.80.Pj, 34.10.+X 



Precise knowledge of atom-atom interactions is of sig- 
nificant importance in the creation and manipulation of 
ultracold atoms or molecules as well as in their applica- 
tion to precision measurements. Investigations on light- 
induced atom losses in magneto-optical traps (MOTs), 
where collision processes are altered by the presence of 
near resonant photons have triggered an emerging 
field of cold collisions with laser cooled atoms 0. The 
stability and dynamics of the Bose-Einstein condensates 
(BECs) are governed by the ultracold coUisional proper- 
ties 

II II nil- Their tunability via magnetic or optical 
method is utilized to produce ultracold molecules as well 
as molecular BECs 0; S IQ]- From a metrological 
point of view, atomic collisions introduce unwanted fre- 
quency shifts that critically limit the performance of the 
state-of-the-art atomic fountain clocks ^-^id optical 
clocks 113. 

A two-body collision problem is most simplified for ul- 
tracold atoms. The ground-state wave function for the 
angular momentum I = (s-wave) state can be written 
as R{r) = Ug{r)/r with Ug{r) the solution of an ordi- 
nary one-dimensional (ID) Schrodinger equation and r 
the interatomic separation. The interaction is character- 
ized by a single parameter, the s-wave scattering length 
a, which is associated with the last node of the scattering 
wave function %(r) 0. One of the most versatile ex- 
perimental techniques for probing the wave function, es- 
pecially its last node, is the photoassociation (PA) spec- 
troscopy Q that measures the intensity modulation of 
the PA lines [isl IT^ . Since the demonstration of the 
technique in laser-cooled Na J8J and Rb 19J atoms, there 
have been extensive theoretical as well as experimental 
studies in alkali-metal systems 0. However, analysis of 
these PA lines are not easy because the hyperfine struc- 
ture in their ground state as well as the excited state gives 
rise to complicated molecular potentials often referred to 
as hyperfine spaghetti 20]. In contrast, alkaline-earth- 
metal system offers the simplest energy structures, i.e., 
the ^Sq ground state and the ^Pi excited state, which 
allows straightforward interpretation of PA spectra. Re- 
cently a PA line profile was obtained for Ca, which de- 
termined possible range of the scattering length _21^]. A 
predissociation process was revealed thr oug h the broad- 
ening of the PA line profiles in Yb atoms W^, ■ PA spectra 



near the dissociation limit was investigated for Sr, which 
determined the ^Pi excited state lifetime 23]. 

In this Rapid Communication we report on the deter- 
mination of the last node of the scattering wave function 
for the 5s^ ^Sq state of ®*Sr atoms by reconstructing the 
wave function through PA spectroscopy. The PA spectra 
were investigated for detunings up to A = —600 GHz be- 
low the 5s^ ^So+SsSp ^Pi atomic asymptote. Thanks to a 
simple electronic structure of alkaline-earth-metal atoms, 
the intensity modulation of the PA spectra straightfor- 
wardly inferred the squared scattering wave function for 
the interatomic separation of 3.2 — 11 nm. The last 
node of the scattering wave function was found to be 
To = 3.78(18) nm, which suggested a possible range for 
the scattering length a. The PA resonances were also 
used to derive the most precise lifetime of 5.263(4) ns 
yet reported for the 5s5p ^Pi state, which will provide a 
stringent test for a many-body calculation of the atomic 
structures [Hill El- 

We consider an atom pair in the ^Sq ground state 
with relatively large interatomic distance r, where the 
resonant dipole-dipole (RDD) interaction dominates the 
atomic interaction. Figure ^ shows schematic energy lev- 
els. By shining a PA laser with frequency which is 
negatively detuned {A = i'l — vq < 0) in respect to the 
^5*0 Pi atomic transition, the atom pair is excited to 
the molecular state bound by the RDD interaction. 
The molecular potential is given by 



'2^ 



(1) 



for r < Ao where the retardation effects |23 are neg- 
ligible. Here D ~ hv^ is the energy of the ^S'o + ^Pi 
atomic asymptote, and h — hjl-K is Planck's constant. 
The coupling constant C3 for the ^11+ molecular state is 
expressed in terms of the reduced wavelength Aq = Ao/27r 
of the corresponding atomic transition and the radiative 
lifetime tq of the ^P\ state. In the following discussion, we 
focus on atomic separations where the van der Waals in- 
teraction energy Vgir) is much less than the excited state 
width of ~ fi7o, where 70 = I/tq is the atomic lincwidth. 
The resonant excitation condition for the molecular po- 
tential is, therefore, given by 



/lA = 14 (r) - Vg{r) - huo ~ -C^/r^ 



(2) 



2 



Energy 



PA Signal 




x atomic separation r 



FIG. 1: Relevant Sr2 molecular energy levels: The 
ground state and the attractive molecular states near 

the 5s^ + 5s5p ^Pi atomic asymptote are shown by dashed 
lines. At the Condon radius rc, the PA laser excites an atom 
pair in the ground state to the fth vibrational levels of the 
^Sj state with a binding energy of /lA. The PA signal in- 
tensity, as schematically depicted in the left panel, probes the 
squared amplitude of the scattering wave function Ug(r) (solid 
line) to determine its nodal point of tq. 



In the spontaneous dissociation of the molecular states, 
the atom pair will acquire enough kinetic energy from 
the molecular potential through radiative decay or state 
changing processes llj to be ejected from a shallow trap- 
ping potential. These atom-loss mechanisms are used to 
detect the PA process. 

The PA spectroscopy was done for a few microkelvins 
cold ^^Sr atoms trapped in a ID optical lattice jl^, which 
was loaded from a narrow-line MOT operated on the 
^5*0 - ^Pi transition at A = 689 nm with a linewidth 
of 7.5 kHz. The ID optical lattice consists of a standing 
wave of a far-off-resonant trap laser with a wavelength 
of Al « 810 nm. Its power of 250 mW was focused to 
the beam radius of 25 ^m, which gave radial and ax- 
ial trap frequencies of 350 Hz and 40 kHz, respectively. 
With this trapping laser parameters, the population of 
the ^Pi state, which is the most strongly coupled to the 
ground state by the trap laser, is approximately 10"^ 
and was neglected. In addition, the induced light shift 
of 5U /h Ri 150 kHz and the Doppler shifts (motional 
sidebands) were much less than the molecular linewidth 
of •jq/tt « 60 MHz. The typical number of atoms A^o 
loaded into the optical lattice was 10^, corresponding to 
a peak density of hq « 10"'^^ cm~^. This high atomic den- 
sity offered us a clear PA signal, as the two-body atom 
loss rate (3no could be made larger than the linear loss 
rate P due to background gas collisions. The ultracold 
temperature of atoms suppresses the higher order partial 
waves and a single rotational level (J = 1) is allowed for 
the state. 

The PA laser at Aq = 461 nm was generated by fre- 
quency doubling of a frequency-stabilized titanium sap- 



phire laser (Coherent, MBRllO), in a 7-mm-long waveg- 
uide periodically poled lithium niobate (PPLN) crystal. 
The mode-matched infrared power of 200 mW was con- 
verted to the blue output power of 70 mW. The PA 
beam with its radius of 120 fj,m was superimposed on 
the ID lattice beam. Its intensity /pa was typically 10 
W/cm^, which is much larger than the saturation inten- 
sity /sat = 42 mW/cm'^ of the associated atomic transi- 
tion. Such a rather intense PA laser was needed to excite 
deeply bound vibrational levels with a binding energy of 
hA, since the Franck-Condon factor decreases as |A|~-'^/^ 
[2^. Because of this A-dependent PA rate, we set a du- 
ration rpA of the PA laser in the range of 5 — 600 ms so as 
to obtain typical atom losses of 50% during a frequency 
scan. The PA laser frequency was monitored by a wave 
meter (Burleigh, WA-1500), which was calibrated by the 
absorption spectrum of Te2 30] . One cycle of the mea- 
surement took 1 s, and the sweep rate of the PA laser 
frequency was 40 MHz/s, which was continuously swept 
over 48 GHz in a single scan. 

The inset of Fig. |21 shows a typical PA atom-loss 
spectrum at A w -107 GHz with /pa = 10 W/cm^ 
and TpA ~ 15 ms. A loss fraction of K(A,rpA) = 
[Nq — iV(rpA)]/A'^o w 0.5 was obtained, where A^o and 
N{tpa) are the number of atoms before and after the 
PA process, respectively, measured by laser induced flu- 
orescence on the ^Sq — ^Pi cyclic transition. This PA 
linewidth 7pa of a few hundred MHz was explained in 
terms of the power broadening, as it was decreased down 
to (84 ± 12) MHz, which is roughly a molecular linewidth 
of 7o/7r w 60 MHz, by reducing the PA laser intensity 
to 40 mW/cm^. By increasing the detuning |A| of the 
PA laser, the PA signal k decreased faster than |A|^^/^ 
and was finally buried in the shot-to-shot atom number 
fluctuation of «10% at A « -250 GHz. We attributed 
this disappearance of the PA signal to the node of the 
ground-state wave function. The vibrational lines reap- 
peared at A « —400 GHz and the PA signal k increased 
with further detuning. The upper panel of Fig. |5| shows 
the change of the PA Hnewidth for /pa = 10 W/cm^. A 
line broadening started from A « —400 GHz, which may 
be attributed to predissociation plj . 

We first analyze the observed PA resonances to confirm 
the validity of assuming the RDD interactions of Eq. 
in our scanning range of —600 GHz < A < —12 GHz. 
For this interaction potential close to the dissociation 
limit, the vibrational energy E{v) can be given by the 
semiclassical formula Isil . 



Eiv) = D - Xo{v - vdY- 



with 
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Here v denotes the vibrational quantum number counted 
from the dissociation limit (the highest bound state is 
w = l),0<'!Ji3<lis the effective quantum number of 
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FIG. 2: PA resonances |Apa(u)| as a function of the vibra- 
tional quantum numbers v plotted on a log-log scale are shown 
in the lower panel. The solid line is the fit |A(w)| by the semi- 
classical formula given in Eq. The inset shows a typical 
PA line; atom loss fraction function of detuning A. 

The residuals of the fit, S{v) = Apa(i') — A(«), are plotted in 
the middle panel. The upper panel shows the PA linewidth 

7PA. 



the dissociation limit, and /i the reduced mass. Figure|21 
summarizes the observed PA lines |Apa(w)| as a function 
of the vibrational quantum number v. These resonances 
were fitted by 



A{v)^~Xoiv--VD)Vh, 



(5) 



with Xq and vd being the fit parameters. The mid- 
dle panel shows residuals 6{v) — Apa(w) — A(w) of the 
fit. The residuals S{v) were comparable to the measure- 
ment uncertainty of the PA laser frequency of « 300 
MHz. Any systematic deviation of the fit was not found, 
which confirmed the validity of Eq. Q that neglected 
the higher order dispersion terms such as Cg and Cg. 
The Xq obtained by the fitting determined the coupling 
coefficient C3 = (18.360 ± 0.013) (atomic unit). This in- 
ferred the atomic radiative lifetime of the ^Pi state to 
be To = (5.263 ± 0.004) ns or the atomic linewidth of 
7o/27r = (30.24 ± 0.02) MHz, which is one order of mag- 
nitude improvement in uncertainty compared with a pre- 
vious measurement p^ . 

The intensity modulation of the PA signal gives infor- 
mation on the amplitude of the ground-state scattering 
wave function Ug{r), as the PA signal intensity is asso- 
ciated with the free-bound Franck-Condon (PC) factor 
F(A). Following the reflection approximation the 
FC factor is calculated by the squared ground-state scat- 
tering wave function evaluated at the Condon radius re. 



F{A) 



dE{v) 1 
dv dc 



\ug{rc)Y 



(6) 



where dc = \{d/dr) [Ve{r) - Vg{r)]\^^^^ is the slope of 
the differential potential at rc- Neglecting the ground- 
state van der Waals potential, which is as small as 



Vg/hK, —200 kHz even at the Condon radius correspond- 
ing to A = —600 GHz, the Condon radius for the uth vi- 



brational state is calculated as rc{v) = yC^/\h\K(v)^ 
using Eq. 0. Assuming the PA loss rate coefficient /3(A) 
is proportional to the FC factor, Eq. lO gives the squared 
ground-state wave function as 



dv 



(7) 



which indicates that the scattering wave function can be 
reconstructed by the PA loss rate sampled at the Condon 
radius. 

The PA loss rate coefficient /3(A) is determined from 
the PA loss fraction k(A, tpa) as 



/3(A) 



r K(A,TpA)/?7(rpA) - 1 



no 



1 - K(A,rpA) 



(8) 



by solving the rate equation for a time-dependent atom 
density n{t) in the presence of collisions, i.e., n — —Tn — 
/3(A)n^. Here ng is the initial atomic density, P = (0.38± 
0.05) sT^ is the linear trap decay rate due to background 
gas collisions, and ri{TpfC) = 1 — e~^'^^* is the linear loss 
fraction. Figure U] shows the reconstructed wave function 
|ug(r)|^ from the PA line intensities using Eq. 0. In this 
way, the disappearance of the PA lines is related to the 
nodal point of the scattering wave function. Assuming 
the ground-state scattering wave function to be linear 
around the node r = rg, i.e., Ug{r) cx (r- ro), the squared 
wave function is given by 



\ug{Tf^A{r~rQf 



(9) 



with A a constant. The experimentally derived squared 
wave function as shown by dots in Fig. |31 was fit by the 
above function, where we weighted the data points by 
their uncertainty. The node of the ground-state wave 
function was determined to be tq = (3.78 ± 0.18) nm. 
Our previous measurements on the rethermalization time 
for ultracold ^^Sr atoms 32] inferred that the scattering 
length is not larger than 5 nm, or the last node to be 
To < 5 nm, which suggests that the observed node is the 
last one. 

The last node position rp can be straightforwardly as- 
sociated with the s-wave scattering length a, especially 
when the node position is larger than a quantum thresh- 
old radius of rg as in the case of Na [31 . For 
®*Sr, however, the observed node position was slightly 
smaller than this characteristic distance of rq = 4.0 nm, 
therefore, a naive usage of the semiclassical formulas 
for calculating the scattering length from the node may 
be problematic. For example, the semiclassical formula 
Eq. (11) in Ref. J7J inferred a positive scattering length 
of a = (1.3 ± 0.5) nm, while a quantum formula Eq. (14) 
in the same reference gives a negative scattering length 
of a = -2tl nm |3|. 

In conclusion, we investigated the PA spectra of ul- 
tracold ^*Sr for detunings up to 600 GHz to the red of 
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FIG. 3: The squared ground-state wave function |u9(r)p as a 
function of the interatomic distance r. Filled squares indicate 
reconstructed wave function from the PA signals. The inset 
is the enlarged view around the node, where the data points 
were used to be fit by = A{r — r\)f (solid line) to 

determine the node to be ro = (3.78 ± O.f 8) nm. 



the dissociation limit. The observed PA lines determined 
the most precise natural lifetime of the 5s5p ^Pi state of 
®*Sr yet reported. The disappearance of the PA spectra 
gave us information on the nodal point of the ground- 
state wave function. The node, unfortunately, did not 
lead to a unique determination of the s-wave scattering 
length in the framework of the semiclassical treatment 
and urged careful quantum treatment. However, small 
positive or negative scattering length suggested by this 
work indicates the difficulty in attaining BEC through 
evaporative cooling in ®^Sr. 

Recently, we became aware of Ref. j35| which measured 
the last node position of the ground-state ®*Sr collision 
wave function. 
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